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The IKKT or IIB matrix model has been postulated to be a non perturbative definition of super- 
string theory. It has the attractive feature that spacetime is dynamically generated, which makes 
possible the scenario of dynamical compactification of extra dimensions, which in the Euclidean 
model manifests by spontaneously breaking the SO(10) rotational invariance (SSB). In this work 
we study using Monte Carlo simulations the 6 dimensional version of the Euclidean IIB matrix 
model. Simulations are found to be plagued by a strong complex action problem and the factor- 
ization method is used for effective sampling and computing expectation values of the extent of 
spacetime in various dimensions. Our results are consistent with calculations using the Gaussian 
Expansion method which predict SSB to SO(3) symmetric vacua, a finite universal extent of the 
compactified dimensions and finite spacetime volume. 
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1. Introduction 

Superstring theory is a natural candidate of a unified theory of all interactions, including grav- 
ity. It lacks a non-perturbative definition, which would allow us to address dynamically questions 
such as the preferred vacuum where our universe sits and determine its properties such as classi- 
cal spacetime dimensionality. The large N limit of the 10 dimensional IKKT or IIB matrix model 
has been proposed to provide such a definition Using dualities, which suggest that the five 
types of superstring theory are actually limits of a unique theory, the IIB matrix model is expected 
to describe the unique underlying theory despite its explicit connection to perturbative type IIB 
superstring theory. 

The IIB matrix model has a series of attractive features: The model has a unique scale, which 
raises the possibility for the theory to choose a unique vacuum. Spacetime and matter content can 
arise dynamically from the distribution of the eigenvalues of the bosonic matrices, which makes 
possible the realization of the scenario of dynamical compactification of extra dimensions. The 
motivation for studying such a scenario in the Euclidean model comes from noticing that lower 
dimensional configurations are stationary points of the fluctuating complex part of the action of the 
model [0]. Using the Gaussian Expansion Method (GEM), one can show that the SO(3) symmetric 
vacuum has the smallest free energy density compared to higher or lower dimensional configura- 
tions [Q]. Dynamical compactification, therefore, occurs by spontaneously breaking the SO(10) 
rotational symmetry of the model (SSB). Moreover, a universal scale r for the small dimensions is 
computed 1 , with the scale R of the large dimensions being such that R d r i0 ~ d = £ 10 , with £ being the 
scale of the size of the symmetric configurations. The former relation is interpreted as a constant 
volume property of the d dimensional (metastable) vacua of the model and leads to the conclusion 
that spacetime in the Euclidean IIB model has finite volume 

Monte Carlo simulations of the Euclidean IIB matrix model, as well as simpler related mod- 
els, can confirm these results from first principles and shine light into the mechanism of dynamical 
compactification of extra dimensions [Q, Q], Such simulations are hindered by a severe complex 
action problem and are quite difficult. Some results, however, have been obtained by studying 
related lower dimensional toy models using the so called factorization method. This method, pro- 
posed originally in [||] has been tested also in Random Matrix theory and finite density QCD [Q]. 
It attempts to sample the most important configurations that contribute to the partition function, 
first by numerically solving equations that compute the maxima of the distribution functions of 
appropriately chosen observables that are strongly correlated to the fluctuating phase ||] and then 
by simulating a constrained system in the region of the solutions. This makes importance sampling 
possible by also taking into account the suppression of configurations by the fluctuations of the 
phase, together with the suppression caused by the measure and the real part of the action in the 
partition function. The overlap problem is thus solved and the complex action problem reduced 
due to milder fluctuations of the phase within the sampled configurations and the use of scaling 
properties that permit extrapolations to larger systems. 

In this work, we present Monte Carlo simulation results of a 6 dimensional version of the 
IIB Matrix Model. Contrary to where the oneloop approximation was used, the full model is 
simulated. This is necessary since, as we will show, the short distance non perturbative dynamics 

'Note that the actual length of each dimension is, by definition, r 1 ' 2 and R l l 2 respectively. 
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of the eigenvalues of the matrices play a crucial role in generating the scales of dynamical com- 
pactification. We show that in the absence of the complex part of the action, no SSB of the SO(6) 
rotational invariance occurs. We apply the factorization method in order to compute the eigenval- 
ues X\ > %2 > • • • > ^6 of spacetime's "moment of inertia" tensor T^ v and compare our results with 
those obtained in ^ using the GEM. In particular, we test the GEM predictions of SSB to SO(3) 
symmetric vacua, a universal scale r for compactified dimensions and finite spacetime volume. Our 
results are consistent with these scenaria. 



2. The Model 



We study a 6 dimensional version of the IIB matrix model defined by the partition function 



Z 



J dAdydyeT Sb - s f , (2.1) 



where Sb = — ^tr [A^,A V ] 2 and 5/ = — tr^a^)^ [A^, y/jj] . The model can be formally obtained 
by the dimensional reduction of the Euclidean 6 dimensional ,jV = 1 \J(N) Super Yang Mills the- 
ory to zero dimensions: The N x N matrices A„, fi = 1, ... ,6 are traceless and hermitian and 
transform as 0(6) vectors. The N x N matrices \j/ a , \j/ a are traceless with grassmannian entries 
and transform as Weyl spinors. The model turns out to have JV = 2 supersymmetry, which leads 
to the interpretation of the eigenvalues of defining the points of the 6 dimensional (Euclidean) 
spacetime 2 [|IJ]. Dynamical compactification can manifest if the distribution of these points break 
0(6) rotational invariance spontaneously (SSB). The order parameters of such SSB are the expec- 
tation values (X\ ), . . . , (A^) of the "moment of inertia tensor" of spacetime T^ v = — tiA^A v , where 
the eigenvalues X\ > . . . > are ordered before taking the expectation value. SSB occurs if, in 
the large ,/V limit, some of the expectation values ( X\ },..., (A^) grow "large" and the remaining 
(Arf+i), . . . , (Xe) remain small. The large eigenvalues define the extended dimensions of spacetime 
and we obtain dynamically a d dimensional spacetime. This scenario has been studied using GEM 
in [||] where S0(3) invariant configurations were found to have the minimum free energy density, 



thus dominating the path integral (2.1). The 6 — d small dimensions in the SO(d) vacua turn out to 



have a ^-independent thickness r 1 / 2 , whereas the d large ones have thickness R l l 2 , such that 

R d r 6 - d =£ 6 , (2.2) 

where I pa 0.627 is the square of the extent of the SO(6) symmetric configurations, (Ai) = ... = 
(X<i) = R and (Aj + i) = . . . = (Xq) = r. The SO(6) symmetric configurations dominate in the phase 



quenched model discussed below where SSB does not occur [g]. Eq. ( |2.2| ) expresses the constant 
volume property of the SO(<i) vacua and implies that the volume of spacetime in the 6 dimensional 
Euclidean IIB matrix model infinite [J3j] . 



Monte Carlo simulations are performed by integrating out the fermions in (2.1) first, obtaining 



J dAdet^e- s "= J dAe,- So+ir , (2.3) 



2 Albeit a non classical spacetime with a fuzzy geometry since the dominant configurations of cannot be simul- 
taneously diagonalizable pp. 
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where det..# = Zf[A] = J dydxj/e Sf and So = Sb — Iog|det.y#|. The determinant det^ = 

|det./#|e ,r turns out to be generically complex creating a very strong complex action problem 
in the simulations. We applied the factorization method [Q], which first amounts to considering the 
phase quenched model Zq = J dAe~ So . In this model, we computed the phase quenched expecta- 
tion values (Ai)o, • • . , (X n )o, and defined the normalized eigenvalues X n = X n j (A„}o- Deviation of 
(X n ) from 1, is the result of the suppression of the dominant configurations of Z by the fluctuations 
of the phase F. Due to the strong correlations of Xn with T, we consider the distribution functions 

M 

p n (x) = (8{x-~X n )) = ^pf\x)w n {x), (2.4) 

where 3 pi (x) = (8(x — X n ))o, w n (x) = (e' r )n, x and (-}, ux denotes expectation values within the 
constrained system Z HjX = J dAe~ So 8(x — X n ). 

The value of (A,,) in the large N limit is determined by the minimum of the "free energy" 



^ n (x) = — logp„(x) which, given Eq. (2.4), is the large ,/V limit of a solution to 



aM 0) w " ^ logpi0) {x) = -W^Tx XogWn{x) ■ (2 " 5) 



Using the dominant solution of (2.5) as an estimator of (X n ) for finite ,/V, has the advantages that 



the overlap problem is solved by simulating Z„ v and the complex action problem reduced, since 
<£>n(x) = lim — ^logw„(x) is found to scale for relatively small ,/V and it is possible to extrapo- 
late solutions to larger values of N. Moreover, the numerical errors of the solution do not grow 
exponentially with ,/V [||, ^]. 

3. Simulations and Results 

Monte Carlo simulations are performed on the system Z n y = J dAe.~ So ~ v (^ where V(z) = 

\y{z — §) 2 and 7, are parameters. The rational hybrid Monte Carlo method (RHMC) is used 
in the simulations. We take y large enough, so that p n y (x) is sharply peaked at x p and the re- 
sults are independent of y. We use the estimators x p = (A„)„y, w„(x p ) = (e' r )„ y and fl°\x p ) = 

(X n ) V'((X n ) n y). 

For 7 = we obtain the phase quenched model Zq. We simulate this system and show the 
results for (X„)q, n = 1, ... ,6 in the left plot of Fig. [T[ In the large ,/V limit, these values converge 
to the same value £ 0.627 as predicted by GEM [ph. No SSB of SO(6) occurs in the absence of a 
fluctuating phase T as expected ^j]. 

In order to compute the solution of (Oh in the large N limit, it is important to use the scaling 



properties of the functions w n (x) and /„ (x) for small x and in the large ,/V limit. We find that <& n (jc), 

n = 3,4,5,6 is almost constant for large x. Thus the extended direction becomes decorrelated 

with the phase and there is no need to constrain the large eigenvalues as in [§]. We also observe 
1 



good scaling of —=\ogw n {x) for 12 < N < 24. For solving (2.5), the small x scaling behavior M\ 
N — 



3 The constant C = (e' r )o is not necessary in the calculations described below. 
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Figure 1: (LEFT) The (A„)o of the phase quenched model are plotted against \/N for 12 < N < 48. The 

/ 6 \ 1/6 

horizontal line is I = 0.627. (RIGHT) h\ = ill (k)n,V J for = 32. The horizontal line is I = 0.627. 

^2 log w„(x) ~ — co, n — ci ifI x 7 ~", « = 2, . . . , 6 is important. We compute the coefficients cq m and 
ci >n , as shown in the left plot of Fig. ^[ For fn°\x) the expected small x scaling is —f„°\x) ~ 

-(7 — n) + 25„i 1 - and implies the existence of a hard core potential suppressing the small x 
2 J x 

region, as it was first found in [^j. This is a strictly non perturbative effect which is absent e.g. 
in the one loop approximation [[]]. Around x = 1 a — /« (x) scaling is found to dominate 4 . This 
adds finite size effects to the small-* scaling which can be subtracted by removing the ff(l/N) 
term g n (x)/Nx, where we take g n {x) = a\, n (x — 1) + a2,„(x — l) 2 + #3,n(x — l) 3 around the peak 
x « 1. The computation of g n (x) is done by the fit in the right plot of Fig. where we show 

the result for n = 4. Finally, the left hand side of (2^) is estimated by the function jpfn^ix) — 

^" = \-(l — n) + 28 n \ 1 CX ^ - — ( k^L^ where q n is determined by the fit of the left plot of Fig. 
Nx [2 J x 

[3| In the same plot we also show —<& n (x), and the intersection of the two curves estimates the 
solution of ( p3| ) for n = 4. It is found that (A4} = 0.31(2). The corresponding GEM prediction is 
rjl rj 0.223/0.627 = 0.355 [g]. Preliminary results for = 3,5 and 6 yield similar results which 
are consistent with the GEM finding that r takes a universal, ^-independent value. The constant 
volume property is further studied by computing L 2 , = (nf=i (^4)n,v) The results are shown in 
the right plot of Fig. [!]. For 0.5 < x < 1 h\ « £ « 0.627 as predicted by the GEM. 

Finally we attempt to compare the free energies J^3(x) and ^(x) by numerically computing 
A34 = — ^(x) + ^(x). Their value at the solution of Eq. ( |2.5[ ) compares the free energies of the 
SO(3) and SO(2) vacua respectively. The calculation is done as described in [||] by computing 



A34 = ^ logw 3 (x. v ) -_logw 4 (x,)-y so(2) — — logp(°)(x) ~ — logw 3 (x,) - — logw 4 (x. v ) 

(3.1) 



Note that in the one loop approximation, the — /„ (x) scaling holds /or all values of x |g] . 
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Figure 2: (LEFT) Small x scaling of log W4(x)/N 2 , according to the discussion on page | for N = 24. 
(RIGHT) —f^°\x) scaling around the peakx « 1, according to the discussion on page |forjV = 24. The fit 
is the function g4 (x) described in the text. 




Figure 3: (LEFT) The solution of Eq. (2.5), giving (A4} =0.31(2). The vertical axis and the corresponding 
fit of the data points estimate the LHS of ( 231) in the large N limit. (RIGHT) log w„ (x) /N 2 for N = 24 with 
the corresponding fits to their small x scaling. These curves can be used to compute the difference A34 of 
Eq. (ph. 



at x s « 0.31. The term / — - — \ogp^\x) vanishes in the large ,/V limit. The comparison 

7so(2) N 1 dx 

of the free energies can be read off the right plot of Fig. ^ where the N = 24 results are shown. 
Although we cannot draw a definite conclusion, we see no inconsistency with a dominating SO(3) 
vacuum as predicted by GEM. 

In summary, we simulated from first principles a 6 dimensional version of the IIB matrix 
model. We studied the scenario of dynamical compactification by SSB of the eigenvalues of the 
tensor T^ v and found that our results are consistent with the quantitative predictions of the GEM. 
The importance of the fluctuating phase in inducing SSB was demonstrated first by confirming 
that no SSB occurs in the phase quenched model. A strong complex action problem was found 
to hinder ordinary Monte Carlo simulations and the factorization method was used in order to 
determine and sample the region favored by the competing effects of the fluctuating phase, the real 



6 



SUSY matrix model of dynamical compactification 



Konstantinos N. Anagnostopoulos 



part of the action and the density of states. The distribution function p n (x) of the eigenvalues (X n ) 
were considered and by using the scaling properties of the factors /„ (x) and w n (x) we computed 
the scale r of the compactified dimensions. The results were found to be consistent with the unique 

^-independent numerical value computed using the GEM and with the constant volume property 

R d r 6-d = e 6_ A st 

non perturbative, hard core potential against the collapse of the eigenvalues 
is shown to be generated dynamically which plays a crucial role in obtaining non trivial solutions in 
the large N limit, generating scales r and R that are comparable and giving finite spacetime volume. 

Recent results in however, have shown that it is possible to study the Lorentzian IIB matrix 
model and obtain an expanding, large, 3 dimensional space which arises after a critical time 5 . This 
is a very exciting possibility, which motivates further study of the IIB matrix model as a non 
perturbative definition of superstring theory and as a model for string inspired cosmology. 
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